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The cancer-associated loss of microRNA (miRNA) 
expression leads to a proliferative advantage and 
aggressive behavior through largely unknown mech- 
anisms. Here, we exploit a model system that recapitu- 
lates physiological terminal differentiation and its reversal 
upon oncogene expression to analyze coordinated 
mRNA/miRNA responses. The cell cycle reentry of myo- 
tubes, forced by the El A oncogene, was associated with 
a pattern of mRNA/ miRNA modulation that was largely 
reciprocal to that induced during the differentiation 
of myoblasts into myotubes. The ElA-induced mRNA 



response was preponderantly Retinoblastoma protein 
(Rb)-dependent. Conversely, the miRNA response was 
mostly Rb-independent and exerted through tissue-spe- 
cific factors and Myc. A subset of these miRNAs (miR-1, 
miR-34, miR-22, miR-365, miR-29, miR-1 45, and Let-7) 
was shown to coordinately target Rb-dependent cell cycle 
and DNA replication mRNAs. Thus, a dual level of 
regulation— transcriptional regulation via Rb-E2F and 
posttranscriptional regulation via miRNAs— confers robust- 
ness to cell cycle control and provides a molecular basis to 
understand the role of miRNA subversion in cancer. 



Introduction 

MicroRNAs (miRNAs) are an evolutionarily conserved class of 
small (18-22 nucleotides) noncoding RNAs that have emerged 
in recent years as pivotal regulators of a variety of cellular pro- 
cesses, including differentiation, growth control, and apopto- 
sis (Bushati and Cohen, 2007; Bartel, 2009). miRNAs were 
originally discovered in worms as developmental regulators 
(i.e., lin-4/let-7; Wightman et al., 1993; Reinhart et al., 2000), 
but they have since been found to have a tissue- or lineage- 
specific pattern of expression in higher organisms as well, under- 
lying their direct involvement in tissue differentiation and 
homeostasis in multicellular organisms (Thomson et al., 2004; 
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Landgraf et al., 2007). miRNAs regulate gene expression at a 
posttranscriptional level, with each miRNA being able to target 
several mRNA species (Bartel, 2009). Thus, these small mol- 
ecules are able to establish complex networks of interactions 
to coordinate cellular responses (Stark et al., 2005; Krol et al., 
2010; Mestdagh et al., 2011; Bissels et al., 2011). Deregulation 
of miRNA expression has been associated with multiple dis- 
eases including cancer (Calin and Croce, 2006; Kloosterman 
and Plasterk, 2006), where a global reduction of miRNAs has 
been often observed as a general trait (Lu et al., 2005; Chang 
et al., 2008; Ozen et al., 2008; Kumar et al., 2009). Although 
loss of miRNAs in cancer (and in particular loss of differentia- 
tion-associated miRNAs [DA-miRNAs]) could be consequential 
to the undifferentiated state of advanced tumors, several mo- 
lecular and genetic lines of evidence suggest that the reduction in 
miRNA levels is functional to cell transformation and associated 
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Figure 1 . E1A-TD mRNA signature (GEX-C2C12). (A) Scheme of the El A-TD model. The bar graph shows the quantification of cell-cycle reentry induced 
by El A or the YH47 mutant by BrdU incorporation assay in two independent experiments (error bars indicate mean and standard deviations). (B) The heat 
map shows the 1 ,462 El A-regulated genes at each time point. Gene expression changes during differentiation ("Myotubes") and after the expression of the 
E 1 A YH47/dl928 mutant ("YH47 24 h") are also shown. The ratio between the YH47 mutant and wt-E 1 A regulation was used to distinguish Rb genes (Rb, 
in orange) and YH47 genes (YH47, in green). (C) Functional analysis of the El A-TD mRNA signature. Bars represent the Min-Max p-value ( — logio scale) 
for the most significant (P < 0.001 ) functional classes calculated over all El A genes, or over up- and down-regulated gene classes separately. (D) A series 
of pie charts summarize the E 1 A-TD mRNA signature. The central chart illustrates the transcriptional class of the signature according to the YH47/E1 A ratio. 
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with an aggressive disease phenotype (Kumar et al., 2007; 
Merritt et al., 2008; Martello et al., 2010). To date, however, the 
phenotypic advantages of reduced miRNA expression to tumor 
cell proliferation and progression remain poorly characterized. 

Given the tight relationship between proliferation and dif- 
ferentiation in physiology and in cancer (Croce, 2008; Hanahan 
and Weinberg, 2011), we hypothesized that DA-miRNAs can 
inhibit proliferation, thus cooperating with known regulatory 
gene expression pathways that, under physiological condi- 
tions, keep cell proliferation tightly reined in. To investigate the 
relationship between cancer-promoting proliferative pathways 
and DA-miRNAs, we used a cell model consisting of commit- 
ted progenitors (myoblasts) induced to terminally differentiate 
into myotubes (terminally differentiated [TD] cells) and then 
forced to reenter the cell cycle by the El A adenoviral oncogene 
(Crescenzi et al., 1995). The adenoviral E1A protein behaves as 
a potent oncogene that can revert the tightly controlled state of 
TD cells by inducing proliferation and dedifferentiation path- 
ways, thus mimicking common alterations in human cancer 
(Nicassio et al., 2005; Bianchi et al., 2007). We used this model 
to explore how mRNAs and miRNA components are regulated 
and how they are integrated to effect specific phenotypes. 

Our results indicate that, during terminal differentiation, 
the expression of cell cycle and DNA replication genes is 
repressed by the Retinoblastoma protein (Rb) tumor suppressor 
pathway. This repression is released by the El A oncogene during 
cell cycle reentry. Additionally, we show that a significant frac- 
tion of Rb-dependent genes are also regulated posttranscription- 
ally in an Rb-independent manner, via a subset of DA-miRNAs. 
This dual level of regulation of Rb-dependent messenger RNAs 
(directly via the Rb-E2F axis and indirectly via Rb-independent 
miRNAs) reinforces the postmitotic block that occurs during 
terminal differentiation, and provides a safety mechanism to 
protect against oncogene-induced proliferation. 

Results 

El A induces reexpression of proliferative 
genes through its interaction with Rb 

Proliferating C2C12 myoblasts were induced to terminally dif- 
ferentiate (TD) into myotubes and were then infected with a re- 
combinant adenovirus expressing E1A (Ad-EIA), which induces 
cell cycle reentry and dedifferentiation (Fig. 1 A). We analyzed 
the gene expression program induced during differentiation of 
myoblasts into myotubes (Fig. 1 B, Myotubes), and after E1A 
expression in myotubes (E1A vs. Mock infected myotubes; 
Fig. 1 B, "E1A" ), using the Affymetrix GeneChip 430 platform. 
We found 1,462 genes significantly regulated (corrected P < 0.05) 
by E1A during cell cycle reentry (henceforth the "E1A mRNA 
signature"), of which 1,047 (71.6%) were up-regulated (E1A-UP 
genes) and 415 (28.4%) were down-regulated (E1A-DOWN 
genes). We observed that E1A largely reverts the transcriptional 
program associated with myotube differentiation (911 oppositely 



regulated genes, P < 0.0001 ; Fig. 1 B). Functional classification 
analysis revealed that, although El A-UP genes are significantly 
enriched in proliferative functions (cell cycle and DNA repli- 
cation), E1A-DOWN genes are primarily involved in develop- 
ment, differentiation, and skeletal muscle functions (Fig. 1 C 
and Table S2). 

E1A interferes with the function of the tumor suppressor 
protein Rb and that of the related proteins pl07 and pl30 (col- 
lectively referred to as "pocket proteins"; Lipinski and Jacks, 
1999). To understand the impact of this pathway on the E1A- 
induced transcriptional reprogramming of TD-C2C12 cells, 
we took advantage of the E1A mutant YH47/dl928 (hence- 
forth YH47), which cannot bind to pocket proteins and which 
is unable to force TD cells to reenter the cell cycle (Fig. 1 A; 
Tiainen et al., 1996; Nicassio et al., 2005). By comparing the 
gene expression profiles of TD-C2C12 cells infected with El A 
or YH47, we observed that the regulation of the majority of the 
El A genes was through binding to pocket proteins (9 10 of 1 ,462 
genes, 62%; Fig. 1, B and D; we refer to these genes as "Rb 
genes"). However, although most of the El A-UP genes were 
Rb-regulated (804 of 1,047, 77%, P < 0.0001), the regulation 
of the majority of the E1A-DOWN genes was Rb-independent 
(309 of 415, 76%, P < 0.0001). By unsupervised TRANSFAC 
analysis, we found a significant enrichment for E2F family bind- 
ing sites at gene promoters of El A-UP genes (in particular the 
Rb-dependent ones; Table S3). Instead, the promoters of YH47- 
induced genes (genes induced by E1A and also by YH47, i.e., 
Rb-independent) were devoid of such binding sites. 

We compared our Rb gene list with genes in a published 
genome-wide analysis of Rb-bound promoters identified through 
a chromatin immunoprecipitation (ChIP) sequencing approach 
(Chicas et al., 2010). We observed a significant overlap (P < 
0.0001 ; Fig. 1 D) and a very high correlation (normalized enrich- 
ment score [NES] = 2.69 by gene set enrichment analysis [GSEA] ; 
Fig. 1 E) between the Rb genes (and in particular those up- 
regulated by E1A) and Rb-ChIP bound promoters. Comparable 
results were obtained in a related cellular model, which consists 
of TD primary muscle satellite cells (MSCs) induced to reenter 
the cell cycle by El A expression (Fig. SI). In this latter model, 
we identified Rb-dependent genes through a different indepen- 
dent approach, namely the genetic depletion of a floxed Rb locus 
by Cre recombinase expression (Fig. SI A). 

Collectively, these data establish that the up-regulation of 
proliferative genes by the interaction with the Rb-E2F pathway 
is the major component of the ElA-induced transcriptional 
(mRNA) response during cell-cycle reentry of TD cells. 

El A down-regulates DA-miRNAs 
through an Rb-independent 
transcriptional mechanism 

We analyzed miRNA expression upon terminal differentiation 
of C2C12 cells and during their subsequent ElA-induced cell 
cycle reentry. By measuring the expression of 335 murine miRNAs 



Side charts show the overlap with Rb ChlP-bound promoters by meta-analysis. (E) GSEA was used to correlate the transcriptional classes of El A-TD mRNA 
signature (Rb genes and YH47 genes) with Rb-bound promoters by ChlP-seq (Chicas et al., 2010). Significance (FDR) and NES are shown. 
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Figure 2. E1A-TD miRNA signature (C2C12-miR). (A) A series of pie charts summarize the El A-TD miRNA signature, which consists of differentiation- 
associated (DA-miRNAs, n = 35) and proliferation-associated (PA-miRNAs, n = 21) miRNAs. Each class is subsequently distinguished according to Rb 
pathway responsiveness (Rb, in orange; and YH47, in green). (B) The heat map shows the regulation of miRNAs in the El A-TD signature. Mature miRNAs 
expressed from multiple genomic loci are also shown (i.e., miR-1 33a-l and -1 33a-2) to allow direct comparison with ChlP-seq data. (C) Meta-analysis of 
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by RT-qPCR (see Materials and methods for details), we found 
a total of 56 regulated miRNAs, which were distinguishable 
according to their class of regulation: DA-miRNAs consist 
of miRNAs induced during differentiation and repressed by 
E1A expression (n = 35); proliferation-associated miRNAs 
(PA-miRNAs) instead consist of miRNAs repressed during dif- 
ferentiation and induced upon cell cycle reentry by El A (n = 21; 
Fig. 2, A and B; and Table S4). We refer to this set of 56 miRNAs 
as the E1A-TD miRNA signature. We obtained similar results 
in MSCs (54 DA-miRNAs and 21 PA-miRNAs in MSCs; Fig. S2, 
A and B; and Table S4), and, in fact, we observed a significant 
overlap between the E1A-TD signatures in the two model sys- 
tems (28/35 DA-miRNAs and 14/21 PA-miRNAs, P < 0.0001). 
E1A expression, therefore, has a greater effect on miRNA 
repression than on miRNA induction, something that is particu- 
larly evident in the MSC model, and is reminiscent of the wide 
repression of miRNA expression observed in human cancer 
(Lu et al., 2005; Ozen et al., 2008; Kumar et al., 2009). 

The regulation of miRNA expression during differentiation 
and upon E1A expression occurs at a transcriptional level, as sup- 
ported by several observations. First, those miRNAs, which are 
organized as genomic clusters and thus share the same promoter, 
also shared the same pattern of regulation by E1A (Fig. 2, B and C; 
and Table S4); similarly, the expression of intronic miRNAs, 
which are usually cotranscribed with the exons of a protein- 
coding gene, mostly correlated with the mRNA expression pattern 
of their host genes (21/30, P = 0.0011; Fig. S3 A). Finally, when 
we compared the expression levels of mature miRNAs and their 
respective precursors (pri-premiRNAs) by RT-qPCR, we observed 
concordant regulation, both during differentiation and upon El A 
expression (Fig. S3 B), which indicates that the observed miRNA 
expression patterns were primarily determined by transcription, 
rather than by posttranscriptional levels of regulation. 

We classified miRNAs according to their upstream tran- 
scriptional mechanisms, using the same approach we adopted 
for mRNAs, i.e., by exploiting the E1A-YH47 mutant to dissect 
Rb-dependent and -independent effects. Almost all DA-miRNAs 
were Rb-independent (32/35, P = 0.0008; Fig. 2, A and B). Con- 
versely, the expression of most PA-miRNAs was significantly 
dependent on Rb (17/21, P < 0.0001; Fig. 2, A and B). This latter 
group of miRNAs contained members of the 17-92 cluster and 
the 25-93-106b clusters, which have previously been shown to be 
E2F targets (O'Donnell et al., 2005). Finally, the Rb-independence 
of DA-miRNAs and Rb-dependence of PA-miRNAs were con- 
firmed in the MSC model (34/54 DA-miRNAs, P = 0.0277; and 
19/21 PA-miRNAs, P = 0.0004; Fig. S2, A and B). 

Myc and MyoD regulate the expression 
of DA-miRNAs 

To further characterize the transcriptional regulation of miRNAs 
of the E1A-TD miRNA signature, we examined ChIP sequencing 



(ChlP-seq) databases of transcription factors (TFs) relevant 
to the model system that we used, including Rb (Lipinski and 
Jacks, 1999; Chicas et al., 2010), Myc (Perna et al., 2012), and 
the muscle-specific TF MyoD (Cao et al., 2010). 

We clustered miRNAs of the E1A-TD signature accor- 
ding to the TF occupancy of their promoter regions (see 
Materials and methods for details) and identified five major 
binding groups (Fig. 2 C). Overall, DA-miRNAs were poorly 
bound by Rb, in agreement with their Rb-independence, 
but were frequently bound by MyoD (group a/b) and/or Myc 
(group b/c/d). Conversely, PA-miRNAs, which were mainly 
Rb responsive, displayed a strong binding to Rb (group c') and 
Myc (group c' and d'). Importantly, a similar binding pattern 
for both DA- and PA-miRNAs was found in MSC-derived 
myotubes (Fig. S2, C and D), and almost all the previously 
ChlP-validated TF-promoter interactions were captured by 
our approach (6/6 MyoD. P < 0.0001; 6/8 Myc, P = 0.1790; 
and 2/3 Rb-E2F, P = 0.0240, "ChIP- validated"; Fig. 2 C and 
Fig. S2 C). The most frequently bound TF was Myc, even among 
DA-miRNA promoters (Fig. 2 C and Fig. S2 C). In agreement 
with this observation, the knockdown of Myc in prolifera- 
ting C2C12 (where Myc is highly expressed) caused significant 
up-regulation of several Myc-bound DA-miRNAs (Fig. 2 D), 
which suggests that Myc could be largely responsible for 
DA-miRNA repression. 

□A miRNAs repress myoblast proliferation 
and promote differentiation by antagonizing 
the Rb-ESF pathway 

We examined whether miRNAs play a role in the establishment 
and maintenance of the postmitotic state. We reasoned that 
DA-miRNAs, the main component of our miRNA signature, 
might act as negative regulators of proliferation. We tested this 
hypothesis by overexpressing, in proliferating myoblasts, 10 
DA-miRNAs belonging to the most highly expressed miRNA 
families found in C2C12 myotubes. Of these, seven miRNAs, 
hereafter referred to as "DA-7," caused a reproducible reduc- 
tion in BrdU incorporation (>20% in three independent experi- 
ments, P < 0.05), four being particularly strong in this instance 
(miR-1, miR-365, miR-22, and miR-34, referred to as "DA-4," 
P< 0.01; Fig. 3 A). 

We used gene expression profiling to characterize the tran- 
scripts regulated upon overexpression of each DA-7 miRNA. In 
agreement with previous studies that used a similar approach 
(Baek et al., 2008; Selbach et al., 2008; Guo et al., 2010), each 
miRNA showed a widespread but mild effect on gene expres- 
sion, with thousands of regulated genes being down- and up- 
regulated (Table 1). Although each miRNA affected different 
gene sets, functional classification revealed a parallel effect for 
all DA-7 miRNAs on proliferation and differentiation pathways 
(Fig. 3 B and Table S2). Indeed, by focusing on the El A signature 



the El A-TD miRNA signature using ChlP-seq datasets. Rows show the normalized binding of each TF at the corresponding miRNA promoter. miRNA-TF 
interactions experimentally validated by direct ChIP experiments are also shown ("ChlP-Validated," see Materials and methods for details]. (D] miRNA 
regulation upon Myc knockdown in C2C1 2 myoblasts was evaluated through RT-qPCR. Two different shRNAs (PIKO-Myc-KDl and KD2) were used. Results 
are normalized to myoblasts expressing an shRNA against luciferase (PIKO-shluc). U5A was used as housekeeping gene. A group of MyoD-bound muscle- 
specific miRNAs (miR-1 , -206, and -486) was used as a negative control. The data shown are from a single representative experiment out of two repeats. 
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Figure 3. Identification of biologically relevant DA-miRNAs and characterization of their effect on gene expression. (A) Different DA-miRNAs were trans- 
fected (50 nM) in proliferating myoblasts (C2C12), and the effect on proliferation was evaluated. The bar graph shows the ratio of BrdU incorporating 
cells normalized over control (SCR, scrambled oligo) in three independent experiments. (B) Functional analysis of the DA-7 miRNA-regulated genes. Bars 
represent the Min-Max p-value (— log 10 scale) of each DA-7 gene set for the most significant (P < 0.001) functional classes. (C) Gene expression changes 
of the 1 ,462 El A genes (as from Fig. 1 B) upon DA-7 miRNA overexpression in C2C1 2 myoblasts. (D) The heat map shows enrichment/depletion of Rb 
genes among miRNA down-regulated or up-regulated transcripts. Down-regulated genes are further stratified according to the presence (DOWN_S) or 
absence (DOWN_NS) of seed matches in their 3' UTR region. (E) GSEA was used to correlate Rb ChlP-bound genes (Chicas et al., 2010) with the gene 
expression changes observed upon miRNA overexpression. The enrichment profile (in green) together with the score (NES) and significance (FDR q-val) for 
each gene set analyzed are shown. 
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Table 1 . DA-7 miRNA-regulated genes and miR targets 



miRNA 


Down 


regulated genes 


Up-regulated genes 




miR targets 




n 


Mean ± SD 


n 


Mean ± SD 


n 


Cons. 


Not Cons. 


miR-1 


2 922 


_ a AO + o op 

W 7 ± VJ. Z 7 


3,326 


0.46 ± 0.28 


476 


1 84 


292 


miR-34a 


2,063 


-0.39 ± 0.20 


2,136 


0.40 ± 0.21 


407 


108 


299 


miR-365 


2,569 


-0.46 ± 0.27 


3,142 


0.40 ± 0.21 


384 


66 


318 


miR-22 


2,136 


-0.37 + 0.18 


2,827 


0.41 ± 0.22 


338 


75 


263 


miR-29a 


1,725 


-0.40 ± 0.20 


2,188 


0.37 ±0.18 


413 


234 


179 


miR-1 45 


1,430 


-0.39 ± 0.21 


2,064 


0.39 ± 0.30 


293 


83 


210 


let-7b 


1,991 


-0.39 ± 0.19 


2,022 


0.40 ± 0.22 


355 


199 


156 



The effect of DA-7 miRNAs on gene expression in proliferating myoblasts (C2C1 2) in three independent experiments (see Materials and methods for details) is summa- 
rized. The number (n) of genes classified as miR targets among all the genes down-regulated by each miRNA, classified as "Conserved" (Cons.) or "Not Conserved" 
(Not Cons.) according to TargetScan 5.1 are also shown. 



(as defined in Fig. 1), we observed that a significant number 
of El A-induced genes were down-regulated upon DA-7 miRNA 
overexpression (P < 0.001 for the DA-4 subset and P < 0.01 for 
miR-29a, miR- 145, and let-7b, which compose the remainder of 
DA-7 miRNAs; Fig. 3 C), which suggests that DA-7 miRNAs 
may cause gene expression changes opposite to those induced 
by the El A oncogene. 

Next, we concentrated on Rb genes (910 genes, as de- 
fined in Fig. 1), the major transcriptional component of the E1A 
signature. Rb genes were enriched within the group of tran- 
scripts down-regulated by DA-4 miRNAs and depleted within 
transcripts up-regulated by DA-4 miRNAs (Fig. 3 D). For the 
remaining DA-7 miRNAs (miR-29a, miR- 145, and let-7b), we 
observed only a minor effect. Furthermore, the enrichment was 
independent of the presence of predicted miRNA-mRNA inter- 
actions within the 3' untranslated region (UTR) or the coding 
sequence (CDS; Fig. 3 D and Table S5), which suggests that 
DA-4 miRNAs mainly acted through inhibition of the Rb path- 
way, rather than directly on Rb gene transcripts. Indeed, genes 
repressed by DA-4 miRNAs showed a significant enrichment 
for E2F-family binding sites in their corresponding promoter 
regions (P < 0.01, Table S3) and correlated remarkably well 
with Rb ChlP-bound genes (false discovery rate [FDR] q-value < 
0.001; Fig. 3 E). 

We hypothesized that DA-7 miRNA expression could also 
promote muscle differentiation. Indeed, many muscle genes 
were induced upon the expression of DA-7 miRNAs (Fig. 4), in 
particular by miR-1, miR-34, miR-22, and to a lesser extent by 
miR-29 and miR- 145. Hence, muscle differentiation genes were 
enriched in the transcript set up-regulated by five of the DA-7 
miRNAs (miR-1, -22, -34, -29, and -145, but not miR-365 and 
let-7; Fig. 4 B). 

□A-7 miRNAs jointly repress 
proliferation and accelerate the onset 
of muscle differentiation 

In silico analyses have recently predicted that modules of 
coregulated miRNAs might function coordinately to ensue 
into a specific biological function within the cell (Tsang et al., 
2010). We therefore investigated whether the miRNAs of the 
DA-7 module might cooperate in the regulation of prolifera- 
tion and differentiation. 



A dose-response experiment performed under conditions 
of increasing expression levels of individual DA-7 miRNAs 
revealed that very low concentrations of individual miRNAs 
(1 nM) had no significant effect on DNA replication (Fig. S4, 
A and B). At the same time, however, we noted that combina- 
tions of different miRNAs (at a 1-nM concentration) were able 
to repress cell proliferation, with the DA-7 subset showing a 
significant effect (P < 0.05; Fig. 5 A). 

The simultaneous expression of DA-7 miRNAs repressed 
several proliferation genes and concomitantly induced differ- 
entiation markers either at the RNA or protein level (Fig. 5 B 
and Fig. S4, D and E). Of note, expression changes induced by 
DA-7 were stronger than those observed under conditions of 
overexpression of individual DA-7 miRNAs (Fig. 5 B). 

Furthermore, in a differentiation assay, we observed an ac- 
celerated onset of cell differentiation into multinucleated myotubes 
when cells expressed the full panel of DA-7 miRNAs at low doses 
(1 nM). We did not observe this effect in cells overexpressing indi- 
vidual miRNAs related to muscle differentiation (miR- 1 , Chen 
et al., 2006; and miR-29a, Wang et al., 2008; Fig. 5, C and D). 

These data suggest that, taken as a single unit, the 
DA-7 miRNAs are able to induce cell cycle arrest and pro- 
mote differentiation. 

□A-7 miRNAs directly and coordinately 
target genes involved in the "G1/S 
checkpoint" and DNA replication 

We aimed to clarify the network of miRNA-mRNA inter- 
actions through which DA-7 miRNAs (and in particular DA-4) 
antagonize proliferation and the Rb-E2F pathway. We devel- 
oped a strategy (shown in Fig. 6) to identify direct targets of 
DA-7 miRNAs from our gene expression profiling, based on 
the evidence that mammalian miRNAs act predominantly to 
decrease target mRNA levels (Guo et al., 2010). In brief, we 
selected as candidate targets the mRNAs that (a) contain in 
their 3' UTR at least one miRNA-responsive element (MRE; 
8mer, 7mer-lA, or 7mer-m8), (b) were robustly expressed in 
the context analyzed (C2C12), and (c) were down-regulated 
upon forced miRNA expression in myoblasts. We isolated >300 
targets (miR targets) for each miRNA (Table l and listed in 
Table S6), including several already validated by other studies/ 
approaches (Table S6; Sethupathy et al., 2006; Xiao et al., 2009). 
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Figure 4. Modulation of muscle differentia- 
tion genes by DA-7 miRNAs. A scheme summa- 
rizes the steps that were followed to analyze 
the regulation of muscle differentiation genes 
after the expression of DA-7 in proliferating 
C2C12 myoblasts. (A) One-way hierarchical 
cluster was generated by Cluster 3.0. (B) The 
heat map shows enrichment/depletion of mus- 
cle differentiation genes among miRNA down- 
or up-regulated transcripts. 
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Among miR targets, we found that a significant number of 
ElA-induced genes and Rb genes were targeted by miR-1 
(P < 0.05), miR-22 (P < 0.001), and miR-34a (P < 0.001). Fur- 
thermore, some genes upstream of the Rb-E2F pathway could 
be also targeted, frequently by two or more DA-7 miRNAs (i.e., 
cyclinDs, CDK6, or E2Fs; Table S6). This finding suggests that, 
besides individual effects, DA-7 miRNAs might coordinate the 
coregulation of functionally connected genes to control cell 
cycle and DNA replication (the so-called "co-targeting" effect; 
Tsang et al., 2010). Therefore, we searched for pathways that 
were relevant for proliferation, within which we could also iden- 
tify multiple mRNA-miRNA interactions. We used the Ingenu- 
ity systems pathway analysis software to identify significantly 
enriched biological functions within the Rb gene set (shown in 
Fig. 7). The two most enriched functional categories were "Cell 
Cycle" (206 genes, P = 2.60 x 10" 57 ) and "DNA replication, re- 
combination, and repair" (214, P = 3.65 x 10" 33 ). We combined 
these two gene sets (n = 270) to build a network consisting of 
138 proteins (nodes) and their interactions (edges; Fig. 7 C). 
A significant number of nodes (35 out of 138, 27.5%; P < 0.0001 
as determined by a Fisher's exact test) were miR targets (Table 1 
and Table S6), and roughly one third of the nodes were targeted 
by at least two different DA-7 miRNAs (Fig. 7 C). 

Within this network, we focused on two specific ca- 
nonical pathways relevant to the control of proliferation: the 
Gl/S checkpoint (nine nodes, five of which were targeted by 



DA-7 miRNAs) and "Cell cycle control of chromosomal rep- 
lication" (19 nodes, of which 6 were DA-7 miRNA targets). 
A canonical representation of these two pathways shows that 
both "activators" of the Gl/S checkpoint (positive regulators of 
the pathway, such as cyclins and CDKs; Fig. 8 A) and struc- 
tural or catalytic components of the DNA replication machin- 
ery (MCMs, RPA, Pol A, and Pol D; unpublished data) were 
targeted by DA-7 miRNAs, which suggests a synergistic effect 
of the entire DA-7 class in silencing these two pathways. To mea- 
sure this effect quantitatively, we compared the mean number 
of miRNA-mRNA interactions of the DA-7 subset of miRNAs 
with that of other miRNAs. We observed a significant enrich- 
ment of either target genes (a mean of 5.43 targets per miRNA; 
P = 0.0013) or MREs (an average of 8.14 MREs per miRNA, 
P = 0.0014) within the Gl/S checkpoint pathway (Fig. S5, 
A and B), but not within the "Cell cycle control of chromo- 
somal replication" (targets, P = 0.4095; MRE, P = 0.4035). 
Many of these miRNA-mRNA interactions are predicted to 
be conserved in higher eukaryotes (a mean of 4.90 conserved 
targets per miRNA). 

We generated reporter plasmids for several miRNAs by 
cloning portions of the 3' UTRs encompassing their pre- 
dicted MREs downstream of a luciferase reporter gene; this 
tool allowed us to assess miRNA-mRNA interactions di- 
rectly by measuring relative luciferase activity. In 9 out of 16 
predicted interactions (56%), luciferase activity of the reporter 
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Figure 5. DA-7 miRNAs cooperatively inhibit DNA replication and accelerate the onset of terminal differentiation. (A) DA-7 miRNAs were transfected 
at low doses (1 nM] alone or in combination in cycling C2C12 myoblasts, and the effect on proliferation was evaluated using the Delfia cell prolif- 
eration kit. Bars show the mean and standard deviation over control (SCR, scrambled oligo) in two independent experiments. The asterisk marks a 
significant (P < 0.05) p-value by Student's t test. (B) The effect of DA-7 overexpression (1 nM) was analyzed by RT-qPCR (24 h after transfection) on 
several proliferation and differentiation genes. (C) Myotube (C2C12) formation was assessed upon transfection of DA-7 miRNAs at low doses (1 nM) 
in a 3-d differentiation assay. Differentiation was followed by visual inspection by phase-contrast microscopy. Enlarged/hypertrophic muscle fibers are 
indicated with arrowheads. (D) At day 1 , multinucleated myotubes (>2 nuclei) were counted after immunostaining with MyHC antibody. Insets show a 
1.5x magnification. Bars: (C) 20 urn; (D) 100 urn. 
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Figure 6. Identification of miRNA direct targets ("miR-targets"). (A] The scheme summarizes the steps that were followed to identify direct targets for 
DA-7 miRNAs (referred as "miR targets," see Table S6 for a complete list). (B) Seed enrichment analysis at 3' UTR sequences of miR-1 expression pro- 
file by Sylamer (van Dongen et al., 2008). The input list consisted of all the genes of the array (21 ,608), ordered from the most down-regulated to the 
most up-regulated upon miRNA expression. The two most overrepresented words correspond to the canonical miR-1 7mers (7mer-lA and 7mer-m8, as 
indicated at the top of the panel). (C) Seed enrichment analysis was repeated for each miRNA dataset with 3' UTR, CDS, or 5' UTR sequences. The 
maximal enrichment score (— logio P-value) observed for heptamers of each overexpressed miRNA was used for comparisons and shown in the bar graph. 

(D) Genes containing at least one MRE (7mer words) were downloaded from the Human TargetScan database (Release 5.1; mouse predictions) and split 
according to the expression level in myoblasts (LOW/HIGH based on the 50th percentile of the average raw signal). The bar graph shows the number 
of MRE-containing genes (conserved or not conserved, according to TargetScan) that were down-regulated upon each DA-7 miRNA overexpression. 

(E) Robustly expressed genes (HIGH, as for D) were stratified according to the presence of conserved MRE (CS-TG) or nonconserved MRE (NCS-TG). 
The histogram shows the distribution of up-regulated (UP), down-regulated (DOWN), and not regulated (NOT REG) genes for each class. The overall 
distribution of all genes (Background, BG) is also shown as a reference. A significant response (P < 0.05, Fisher's exact test) was observed for genes 
with either conserved or nonconserved sites. 
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Figure 7. Network of miR targets. (A) The scheme summarizes the steps that were followed to generate a network of DA-7-targeted genes involved in 
the control of proliferation and DNA replication. (B) The bar graph shows the most enriched biological functions as assessed through Ingenuity systems 
pathway analysis software. The list of El A_Rb-induced genes (804] was used as input. A red line marks the threshold for significance. (C) The gene 
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was indeed significantly reduced by transfection of the cor- 
responding DA-7 miRNAs into HEK293T cells (Fig. 8 B 
and Table S7). 

DA-7 miRNAs inhibit El A-induced 
transcriptional reprogramming and the 
induction of cell cycle reentry 

Because DA-7 miRNAs and the El A oncogene had opposing 
effects on the control of proliferation, we reasoned that the abla- 
tion of DA-7 miRNAs in myotubes should enhance cell cycle 
reentry by El A. We tested this hypothesis directly, by knock- 
ing down the strongest DA-7 miRNAs (i.e., the DA-4 miRNAs: 
miR-1, miR-34, miR-365, and miR-22) with anti-miRs and 
measuring cell cycle reentry upon El A expression in differenti- 
ated cells. For all tested miRNAs, we scored a significant and 
reproducible (P < 0.05 by Student's t test; n = 3) increase in 
BrdU-positive cells (Fig. 9 A). Conversely, no effect was observed 
by knocking down miR-24, which also scored negatively in 
overexpression experiments (Fig. 3 A and Fig. 9 A). Therefore, 
DA-4 miRNAs at their physiological levels are able to restrict 
unscheduled cell proliferation induced by the E1A oncogene. 

DA-7 and El A also had opposing effects on gene ex- 
pression, and in particular on the Rb-dependent component 
(Fig. 3 C); thus we analyzed whether endogenous DA-7 miRNAs 
could repress their Rb target genes during cell cycle reentry upon 
E1A expression in differentiated cells. At 24 h after infection, 
when we evaluated the gene expression by El A, many DA-7 
miRNAs were still expressed at higher levels compared with cy- 
cling myoblasts (Fig. S5 C), which indicates that they might still 
be active on their targets. We therefore stratified El A genes by the 
number of MREs for DA-7 miRNAs and observed that mRNAs 
induced by El A were significantly less induced if they inter- 
acted with DA-7 miRNAs (P < 0.0001), with a dose-dependent 
effect (the higher the number of MREs, the lower the induction 
of genes; Fig. 9 B, "E1A on C2C12 myotubes"). This effect was 
not dependent on the length of the 3' UTR regions (adjusted P < 
0.001) and was specifically restricted to Rb genes, as YH47 genes 
were all equally regulated by El A regardless of the presence and 
the number of MREs (Fig. 9 B). Importantly, we were also able to 
reproduce these findings in primary myotubes: upon El A expres- 
sion, the induction of Rb-dependent genes with MREs responsive 
to DA-7 miRNAs was dampened down in a seed-dependent fash- 
ion, whereas no significant effect was observed on Rb-independent 
genes (Fig. 9 B, "E1A on MSC myotubes"). 

These results show that DA-7 miRNAs play a role in pre- 
venting cell-cycle reentry by curbing changes in gene expression 
that are induced by the El A oncogene and confirm that there exists 
an antagonism between DA-7 miRNAs and the Rb pathway, even 
when cells express endogenous levels of miRNAs (Fig. 9 C). 

Discussion 

A coordinated mRMA/miRNA expression 
program in the control of proliferation 
and differentiation 

We have uncovered a coordinated mechanism of control of 
proliferation and differentiation, which couples two distinct 



processes: the transcriptional regulation of cell cycle and DNA 
replication genes, which lies under the control of the Rb-E2F 
pathway, and an miRNA pathway associated with differentiation 
(DA-miRNAs), which targets a significant fraction of the same 
genes, thus conferring robustness to the phenotypic output. 

Our results confirm that the key effect of El A at the tran- 
scriptional level is its interference with Rb and its cognate pro- 
teins (pi 07 and pi 30, "pocket-proteins), which are involved 
in both differentiation and cell cycle withdrawal (Lipinski and 
Jacks, 1999; Ferrari et al., 2008). In parallel, the ElA-associated 
miRNA expression program consists of a widespread repres- 
sion of DA-miRNAs in an Rb-independent fashion. 

During muscle differentiation, the induction of many tissue- 
specific miRNAs (so-called "myomirs") lies under the con- 
trol of muscle-specific regulatory factors (MRFs), such as MyoD 
and myogenin (Williams et al., 2009). We observed an Rb- 
independent regulation of myomirs (miR-1, miR-1 33, miR-206, 
and miR-486), which is in agreement with the known ability 
of E1A to interfere with muscle-specific transcription through 
its N-terminal region, independently of its interaction with 
pocket proteins (Tiainen et al., 1996; Puri et al., 1997; Nicassio 
et al, 2005). 

Our work also demonstrates Myc-dependent repression 
of several DA-miRNAs. Myc is down-regulated during muscle 
differentiation (Yeilding et al., 1998), whereas it is rapidly in- 
duced when cells reenter a proliferative state from quiescence 
or terminal differentiation (for reviews see Amati et al., 1998; 
Grandori et al., 2000). We also observed that many Myc-depen- 
dent miRNAs are induced in serum-starved primary fibroblasts, 
which are quiescent, compared with proliferating cells (unpub- 
lished data). Consistently, Myc overexpression in other cellular 
systems has been shown to repress the transcription of the same 
miRNAs (Chang et al., 2008). Of note, Myc itself is a direct 
target of two miRNAs (miR-34 and miR-145) that are under 
the transcriptional control of p53 (He et al., 2007; Sachdeva 
et al., 2009; Cannell et al., 2010), which is required for complete 
muscle differentiation (Porrello et al., 2000). Hence, it is tempt- 
ing to speculate that the coordinated regulation of miRNAs by 
Myc and p53 might be generally involved in the control of cell 
cycle exit/reentry in nonproliferative contexts (i.e., differentia- 
tion, quiescence, and senescence) as well as in cancer (Chang 
et al., 2008; Kota et al., 2009). 

A module of differentiation associated 
miRNAs act as "quiescence keepers" 

We demonstrated that several DA-miRNAs (the DA-7 miRNAs) 
play a role in the control of proliferation. When overexpressed 
in proliferating myoblasts, where they are normally expressed 
at low levels, these DA-7 miRNAs induce cell-cycle arrest; in- 
stead, the inhibition of the expression of four of the seven DA-7 
miRNAs in myotubes, where they are normally abundantly 
expressed, is sufficient to enhance El A-induced DNA replication. 
These miRNAs reinforce the postmitotic barrier against on- 
cogene-induced proliferation. Therefore, we propose the term 
"quiescence keepers" to describe these miRNAs. Besides indi- 
vidual effects, the quiescence-keeper miRNAs also cooperate 
to induce cell cycle arrest and differentiation. The combination 
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Figure 8. DA-7 miRNAs target Gl/S checkpoint genes. (A) Representation of the Gl /S checkpoint canonical pathway. Genes targeted by at least one of 
the DA-7 miRNA are highlighted in green with a color code (from light to dark green) according to the number of miRNA-mRNA interactions. Those DA-7 
miRNAs targeting each node are reported below in the yellow box. (B) Relative firefly luciferase activity of 3' UTR reporter constructs for some genes of the 
Gl/S checkpoint, after transfection with selected DA-7 synthetic miRNA mimics. Error bars represent standard deviations from three independent transfec- 
tions. Asterisks mark significant values (P < 0.05; Student's f test). Two predicted interactions (CDC25A-miR-l 45 and E2F3-let-7b), for which we did not 
observed a decrease in mRNA levels upon DA-7 expression (and, therefore, were not considered as miR targets), were also analyzed (no-TG). 



of DA-7 had a significant biological effect even when miRNAs 
were expressed at very low (paraphysiological) levels, which 
suggests that the joint regulation of different miRNAs has stron- 
ger activity than individual miRNAs, as observed in previous 
studies (Marasa et al., 2009). 



We identified a set of putative direct gene targets ("miR 
targets") for each quiescence-keeper miRNA, and we built 
a network of miRNA-mRNA interactions that reveals how 
miRNA and gene expression pathways intersect to regulate 
cell cycle exit/reentry. We found that miRNAs target several 
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Figure 9. Endogenous miRNAs inhibit oncogene-induced reprogramming. (A) miRNA inhibitors (anti-miRs] were transfected (100 nM) in differentiated 
myotubes (MSC) infected with suboptimal doses (MOI: 50] of Ad-El A virus. The effect on proliferation was evaluated by means of BrdU incorporation. The 
bar graph shows the ratio of the BrdU-incorporating cell normalized over control (SCR, scrambled oligo) in three independent experiments. Anti-miR-24 
was used as negative control. (B] The effect of endogenous DA-7 on El A-induced genes was evaluated in C2C1 2 (left) and MSC (right) myotubes. El A- 
induced genes were distinguished into two transcriptional classes (Rb and YH47 genes, for C2C12; and Rb-dependent and Rb-independent genes, for 
MSC), and each gene set was further stratified according to the miR target list (Tables 1 and S6) and the number of MREs present in each 3' UTR region. 
Also reported are the mean diamonds (in gray) and the p-value (ANOVA). Three asterisks mark significant distribution (P < 0.05) according to Dunnet's 
t test. (C) A scheme summarizes the role of DA-7 miRNAs to antagonize reactivation of Rb genes involved in cell cycle and DNA replication. 
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key nodes of the Rb-E2F network, including cyclins and Cdks, 
which regulate cell cycle checkpoints and Rb phosphory- 
lation (Sherr, 1996), as well as E2Fs and DPI TFs, which 
mediate the transcription of genes necessary for cell cycle pro- 
gression (Trimarchi and Lees, 2002). In addition, downstream 
effectors of the Rb pathway are also targeted by miRNAs, thus 
revealing multiple levels of interaction between miRNAs and 
Rb-dependent genes that contribute to maintain the balance of 
the postmitotic state. 

In silico analyses have recently predicted that miRNAs 
might function coordinately within a given pathway by cotar- 
geting several genes involved in the same biological function 
(Tsang et al., 2010). We have now demonstrated this same con- 
cept experimentally, showing that functionally related miRNAs 
(the quiescence keepers) act by cotargeting genes of the Gl/S 
transition to reinforce the inhibition of cell cycle progression 
and accelerate the onset of differentiation. 

We believe that our approach, which consists of network- 
ing a "group of co-regulated miRNAs" (with a similar phenotype) 
to a "group of co-targeted genes" (involved in the phenotype), 
can be applied to other experimental settings to reveal func- 
tional mechanisms that would be difficult to uncover with a 
"classical" single gene/miRNA approach. 

From differentiation to tumors: 
could quiescence keepers be 
tumor-suppressor miRNAs? 

Several studies suggest that a reduction in miRNA levels is 
functional to cell transformation, although the phenotypic 
advantages afforded to cells by reduced miRNA expression 
remain unclear (Kumar et al., 2009; Martello et al., 2010). 
Our results suggest that one such advantage may be the loss 
of control of quiescence. In fact, by inhibiting the reexpression 
of oncogene-induced genes in a seed-dependent fashion, quies- 
cence-keeper miRNAs work as a safety net against transforma- 
tion. We have shown how a module of endogenous miRNAs 
can counteract the transcriptional alterations induced by an 
oncogene (E1A, in this case), and we suggest that the miRNA 
pathway may not only serve to stabilize the transcriptional out- 
put of a cell ("buffering" effect), but may also work in prevent- 
ing undesirable alterations because of oncogenic disruption of 
normal cell function ("keeper" effect). 

We speculate that other tissue/context-specific quies- 
cence-keeper miRNAs exist, which are likely to contribute to 
tumor suppression through the regulation of proliferation and 
differentiation. Their identification would provide a spectrum of 
potentially valuable tumor markers. In this regard, we note that 
tumors with low levels of expression of miRNAs tend to be as- 
sociated with a poorer prognosis than those with normal miRNA 
levels (Lu et al., 2005; Merritt et al., 2008), which suggests that 
patient stratification on the basis of expression levels of quies- 
cence-keeper miRNAs may be helpful for prognostic purposes. 

A deeper understanding of the mechanisms responsible 
for determining the balance between differentiation and cell 
proliferation may help to develop new therapeutic strategies in 
cancer. Indeed, the reexpression of quiescence keepers in tu- 
mors can induce cell-cycle arrest and differentiation of cancer 



cells (Kumar et al., 2008; Kota et al., 2009; Taulli et al., 2009; 
Ting et al., 2010), and animal tumor models have been recently 
used to reexpress tumor-suppressive quiescence-keeper miRNAs 
(miR-34, let-7, and miR-1) with promising results (Kim et al., 
2011; Liu et al., 2011; Trang et al., 2011; Wu, 2011). Thus, we 
envision that quiescence keepers might hold a key to new strate- 
gies for differentiation therapy in cancer. 

Materials and methods 

Cell culture procedures 

C2C12 and MSC-Rb lo ' <p/l ° xP myoblasts were differentiated in vitro to myo- 
tubes as described previously (Camarda et al., 2004; Nicassio et al., 
2005). In brief, mouse MSCs were induced to differentiate by plating cells 
(150,000/35 mm) on gelatin-coated dishes (Iwaki), cultured in DME sup- 
plemented with 10% FBS, and, 24 h later, switched to 5% horse serum for 
48 h. Mouse C2C12 myoblasts were induced to differentiate by plating 
the cells (350,000/35 mm) in collagen-coated dishes (Iwaki) and culturing 
them in DME supplemented with 2% horse serum for 3 d. To eliminate 
undifferentiated cells, 50 mM cytosine-D-arabinofuranoside was added 
for the first 48 h. 

Myotubes were infected with an adenovirus carrying the 1 2S form 
of E 1 A (dl520 amino acid changes: del 1 40-1 86; named El A), with the 
YH47/dl928 El A mutant (YH47/dl928, amino acid changes: Y47H and 
C 1 24G; named YH47) or with a control adenovirus expressing a deletion 
of essentially the entire E 1 A gene (dl3 1 2 amino acid changes: del 1 -289; 
named MOCK). 

Multiplicity of infection (MOI) was 300 for all the viruses. The depletion 
of the floxed Rb gene in MSC Rbflox/flox myotubes (provided by A. Berns 
and M. Vooijs, The Netherlands Cancer Institute, Amsterdam, Netherlands) 
was achieved by Ad-CRE infection at MOI 1 ,000 (Camarda et al., 2004). 
Correct differentiation of myotubes was verified by immunophenotypic 
analysis using an antibody against a muscle-specific protein (myosin 
heavy chain [MyHC]) and by RT-qPCR of early and late differentiation 
markers (MyoD and MyoHC genes, respectively). Induction of cell cycle 
reentry by El A was assessed by continuous measurement of BrdU (20 uM) 
incorporation over 24 h. MyHC-positive cells were counted to ensure that 
only fully differentiated cells were considered. Antibodies used in immuno- 
fluorescence: BrdU, clone Bu20a (DakoCytomation); and MyHC, rabbit 
polyclonal or mouse monoclonal (MF-20; a gift from G. Cossu, University 
College London, London, England, UK). Secondary antibodies were pur- 
chased from Molecular Probes or Invitrogen. Nuclei were counterstained 
with Hoechst 33258 or DAPI. Antibodies used in Western blotting: pRb, 
clone G3-245 (BD); p-tubulin, clone 2-28-33 (Sigma-Aldrich); Myc, clone 
Y69 (catalogue no. ab32072; Abeam); PCNA (catalogue no. 610664; 
BD); MyoD, C-20 (catalogue no. SC-304; Santa Cruz Biotechnology, Inc.); 
Myogenin, F5D (catalogue no. SC-12732; Santa Cruz Biotechnology, 
Inc.); NP95 (home-made rabbit polyclonal); and Vinculin (catalogue no. 
V9131; Sigma-Aldrich). Detection was performed with HRP-conjugated 
secondary antibodies (Cappel or Bio-Rad Laboratories) and ECL Western 
blotting reagent (GE Healthcare) or a West Dura kit (Thermo Fisher Scien- 
tific). Chemiluminescence of Vinculin, Myc, PCNA, MyoD, NP95, and 
Myog was performed with the ChemiDoc XRS+ (Bio-Rad Laboratories), 
and band densitometry was analyzed with Image Laboratory 3.0.1 
(Bio-Rad Laboratories). 

Microscopy 

Visual inspection of differentiation was performed under a light microscope 
(DM5500B; Leica) using a 5x objective lens (NA 0.15) at RT and using 
Mowiol mounting medium. BrdU and MyHC staining were viewed under a 
fluorescent microscope (DM5500B; Leica) using a 5x/NA 0.15 objective 
lens at RT. Fluorochromes used were Cy3 (for BrdU) and Alexa Fluor 488 
(for MyHC). All images were acquired with a camera (DFC350FX; Leica) 
and LAS-AF image software (Leica). Images were processed with ImageJ, 
and imaged figures were constructed in Illustrator (Adobe). 

miRNA transfection 

For experiments shown in Fig. 3 and Fig. 4, MiScript miRNA mimics 
(QIAGEN) were transfected at 50 nM into proliferating C2C12 cells, with 
HiPerFect reagent (QIAGEN). After 36 h, RNA was collected for gene 
expression analysis and a BrdU (20 pM) incorporation assay (1 h, 30 min 
pulse) was performed. For experiments shown in Fig. 9, miRNA knockdown 
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in MSC-derived myotubes was obtained by transfection of 1 00 nM miRIDIAN 
miRNA inhibitors (Thermo Fisher Scientific). 1 h before transfection, myotubes 
were infected with the dl520 adenovirus at MOI 50. Cell cycle reentry 
was evaluated by measurement of BrdU (20 uM] incorporation at 40 h 
after infection. 

For experiments shown in Fig. 5, C2C12 cells (p28) were trans- 
fected with control (scrambled) or 1 nM miRNA mimics. For the differen- 
tiation assay, cells were plated on gelatin-coated dishes (400,000 cells/ 
35 mm) and switched (day 1, 24 h after transfection) to differentiation 
medium (2% horse serum). Multinucleated (>2 nuclei) myotubes were as- 
sessed either by visual inspection using phase-contrast microscopy or by 
immunostaining with MyHC antibody (MF20) and DAPI. 

Assessment of BrdU incorporation using the Delfia cell proliferation kit 

Proliferating C2C1 2 cells (p28) were plated in 96-well plates (Viewplate, cata- 
logue no. 6005 181; PerkinElmer) at a concentration of 4,000 cells/well. Cell 
were transfected in triplicate with miScript miRNA mimics (QIAGEN) at a final 
concentration ranging from 0.5 nM to 50 nM in a final volume of 200 ul with 
HiPerFect reagent (QIAGEN) and using the reverse transfection protocol ac- 
cording to manufacturer's instructions. When assessing the combinatorial 
effect of DA-7 miRNAs, mimics were diluted together to give a final concen- 
tration of 1 nM each. At 34 h after transfection, BrdU was added to the 
medium (2 h pulse), and incorporation was measured subsequently using 
the Delfia cell proliferation kit (PerkinElmer) according to the manufacturers' 
instructions and using the EnVision Multilabel Reader (PerkinElmer). 

Luciferase reporter assays 

HEK293T cells (280,000) were plated in a 24-well plate. Shortly after 
plating, cells were transfected with 150 ng of the indicated 3' UTR reporter 
construct (pmiR-GLO; Promega) and miRNA mimics (QIAGEN) at a 50-nM 
final concentration using Lipofectamine 2000 (1.5 ul/well; Invitrogen) 
according to the manufacturer's instructions. After 24 h, cells were lysed 
and assayed in triplicate wells for firefly and renilla luciferase activity 
using the Dual-Luciferase Reporter Assay System (Promega) and an EnVision 
Multilabel plate reader. Firefly luciferase activity was normalized to 
renilla luciferase activity for each transfected well as follows: [firefly(miR)/ 
renilla(miR)]/[firefly(SCR)/renilla(SCR)]. Data are representative of two 
or three independent experiments performed on different days (see 
also Table S7). 

Gene expression analysis 

Gene expression of El A-infected C2C1 2-myotubes was performed through 
the Affymetrix GeneChip Mouse genome 430 2.0 platform. Raw data 
together with detailed description of the normalization procedure is avail- 
able in the GEO database (GSE28457). Significantly regulated genes 
were selected in two independent biological replicates, each performed 
in two technical replicates. The criteria for inclusion in the list of "El A regu- 
lated genes" were: >0.25 log 2 ratio for up-regulated genes and <— 0.25 
log2 ratio for down-regulated genes, and P < 0.05 by Welch's f test with 
Benjamini and Hochberg multiple testing correction. A ratio between the 
regulation of the YH47 mutant and wt-El A was used to derive two classes 
of regulation for El A genes: Rb genes (YH47/E1A < 0.40) and YH47 
genes (YH47/E1A > 0.40). Gene expression of El A-infected MSC-Rb lo,<p/loxP 
myotubes was performed as described previously (Camarda et al., 2004). 
In brief, total mRNA was converted into cDNA and then biotin-labeled 
cRNA using the Superscript II (Invitrogen) and ENZO (Affymetrix) kits. cRNAs 
were hybridized to MGU74Av2 oligonucleotide microarrays (Affymetrix) 
according to manufacturer's instructions. Data analysis was performed 
using Affymetrix Microarray Suite v.5.0 and Excel software (Microsoft). 
E 1 A-regulated genes (>0.25 log 2 ratio for up-regulated genes and <— 0.25 
log2 ratio for down-regulated genes, P < 0.05 by Welch's f test) were clas- 
sified in two classes according to the ratio between Rb~^~ MSC myotubes 
and El A-infected MSC myotubes: Rb-dependent genes (Rb/EIA > 0.40) 
and Rb-independent genes (Rb/El A < 0.40). 

Gene expression profiling upon miRNA overexpression was per- 
formed by Affymetrix Mouse Gene 1 .0 ST arrays. Three biological samples 
were analyzed for each condition. Regulated genes were selected in 
matched experiments as follows: >0.2 mean log 2 ratio for up-regulated 
genes and <— 0.2 mean log 2 ratio for down-regulated genes, P < 0.05 by 
Welch's f test with Benjamini and Hochberg multiple testing correction. 
Data are available in the GEO omnibus database (GSE26764). 

miRNA expression analysis 

The E1A-TD miRNA signature was obtained with a TaqMan Low Density 
Array microRNA Signature Panel V2.0 (Applied Biosystems). Total RNA from 
two independent biological replicas (either for C2C12 or MSC myotubes) 



was used for miRNA profiling. Raw data (Ct values) were exported in 
Excel (Microsoft) for filtering and normalization. Low (mean Ct > 30.01) and 
not expressed miRNAs (i.e., "Not Amplified") were excluded from the 
analysis. Then, median expression of housekeeping genes (U6b, Snol35, 
and Sno202) present on the array was used for data normalization. Regu- 
lated miRNAs were selected based on the following criteria: (a) DA-miRNAs, 
induced upon differentiation (>0.25 log 2 ratio between myotubes and 
myoblasts) and repressed by El A expression (<— 0.25 log 2 ratio between 
El A vs. mock infected myotubes), and (b) PA-miRNAs, repressed upon dif- 
ferentiation (<— 0.25 log 2 ratio between El A vs. mock infected myotubes 
and myoblasts) and induced by El A expression (>0.25 log 2 ratio El A 
infected [El A] vs. mock-infected myotubes). Classification of miRNAs as 
Rb-dependent or Rb-independent was performed as for mRNAs. Several 
miRNAs were analyzed as precursors (pri- and premiRNAs) or mature 
miRNAs by means of RT-qPCR performed with the miScript System (QIAGEN; 
shown in Fig. S3 B). Reactions for mature miRNAs were performed using 
5 ng cDNA for each reaction and U6b/U5a as housekeeping genes. 
miRNA precursors were measured by using 20 ng cDNA for each reaction 
and 1 8s as housekeeping gene. 

GSEA 

GSEA analysis was performed using the preranked gene list (based on 
log 2 expression ratio) of the El A regulated genes (i.e., the 1,462 genes, 
Fig. 1 E) or miRNA-regulated genes (Fig. 3 E). Those genes bound by Rb 
both in quiescent and senescent cells according to ChlP-seq experiments 
(Chicas et al., 2010) were used as gene sets to calculate the enrichment 
score (ES). P-values were calculated by performing 1,000 random permu- 
tations of genes labels to create ES null distribution. 

TF analysis (TRANSFAC) 

E 1 A-TD mRNA signatures from C2C1 2 or MSC myotubes or miR-regulated 
genes (C2C12) were analyzed for TF binding sites, as defined in TRANS- 
FAC version 7.4. Specifically, we used the "compute overlap" tool at 
the Molecular Signatures Database (MSigDB; Subramanian et al., 2005). 
Only significantly enriched gene sets (P < 0.01) matching known TFs 
were considered. 

ChlP-seq meta-analysis of El A genes 

High-confidence Rb-binding regions were previously assessed through 
ChlP-seq in proliferating, quiescent, and senescent human IMR90 fibro- 
blasts (Chicas et al., 2010). Binding regions were mapped to the human 
genome assembly (NCBI36/hg 1 8) by associating each Ref_seq to a pro- 
moter region (considering the 1 ,000 bp upstream or downstream of the 
transcription start site). In particular, we considered those bound both in 
quiescent and senescent cells as Rb-bound genes. Finally, Human Entrez_ 
IDs of these genes were converted into Mouse EntrezJDs to classify El A 
genes or miR-regulated genes as Rb "ChlP-bound" or "Not-bound." 

Genome-wide binding of MyoD, Myc, and Rb to miRNA promoters 

To identify TFs associated with promoters of the El A-TD miRNA signature, 
a three-step procedure was followed. 

First, we retrieved MyoD-, Myc-, and Rb-binding regions, as assessed 
by Chip-seq experiments available from the literature: MyoD (mouse, 
C2C12; Cao et al., 2010), Myc (mouse, 3T9-MEF; Perna et al., 2012), 
and Rb (human, IMR90; Chicas et al., 2010). Only high-confidence bind- 
ing regions, as determined by each original contribution, were used. Then, 
we retrieved miRNA transcription start sites (TSS) along with their TSS 
scores from the work of Marson et al. (2008). TSS scores below zero were 
excluded from subsequent analyses. Of note, we re-annotated the genomic 
coordinates of some well-characterized muscle-specific miRNA promoters 
(miR-1-1, miR-1-2, miR-133a-l, miR-133a-2, miR-206, and miR-486; Rao 
et al., 2006; Small et al., 2010), as data produced in embryonic stem 
cells were likely to be inaccurate in the case of these highly tissue-specific 
miRNAs. All the genomic coordinates were converted to the last genome 
assembly (NCBI Build 36. 1 /hg 1 8 and NCBI Build 37/mm9) using the 
liftOver tool of the UCSC genome browser. Overall, 550 promoters were 
available for human miRNAs and 393 for mouse miRNAs. Because the 
genomic intervals corresponding to each TSS varied from 200 bp to several 
kilobases, a region of 2,000 bp centered on the middle of the TSS region 
was selected ("miR promoter"). Finally, genomic intervals corresponding 
to TF binding and to miR promoters were joined, using mouse promoters 
for MyoD and Myc and human promoters for Rb. Overlapping regions 
(at least 1 bp) were considered as "TF-bound miR promoters." When multi- 
ple peaks were found on a given promoter, only the highest was considered. 
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The height of the highest peak was finally normalized in a 1-100 scale 
to allow their representation using the same scale for different ChlP-seq 
experiments and reported in a log 2 scale (Fig. 2 C and Fig. S2 C). 

ChIP validated miRNA-TF interactions 

Several miRNA-TF interactions shown in Fig. 2 C and Fig. S2 C were 
experimentally validated by direct ChIP experiments ("ChIP validated" 
columns) according to the following references: MyoD (Cao et al., 2006; 
Rao et al., 2006; Small etal., 2010), MYC (O'Donnell etal., 2005; Chang 
et al., 2008; Gao et al., 2009), and E2F1 (Woods et al., 2007; Bueno 
etal., 2010]. 

MYC knockdown in proliferating C2C12 

Experiments were conducted using a pLKO.l lentiviral vector expressing 
two different hairpins complementary to mouse Myc or, as a control, a 
hairpin complementary to luciferase. These vectors were a gift of B. Amati 
(Istituto Europeo di Oncologia, Milan, Italy). The sequences used for MYC 
RNAi were: Myc KD1, 5 '-CTGGAGATGATGACCGAGTTA-3 ' ; and Myc 
KD2, 5 '-CCTGAAGC AGATC AGCAAC AA-3 ' . After four rounds of infec- 
tion, cells were selected with 2 ug/ml puromycin. Total RNA was collected 
after 3 d of selection, and Myc knockdown was evaluated by RT-qPCR 
using a QuantiTect Primer Assay (QT00096194, QIAGEN). 

Seed enrichment analysis by Sylamer 

All the Ref_seq expressed in C2C1 2 cells from the Affymetrix Mouse Gene 
1 .0 ST arrays (21608) were ordered from the most down-regulated to the 
most up-regulated upon overexpression of each miRNA in C2C12 myo- 
blasts. This list was uploaded in Sylamer (van Dongen et al., 2008) to 
compute word enrichment. The analysis was performed with sequences of 
the 3' UTR, which usually contains the MREs, as well as with the CDSs or 
the 5' UTRs for the mouse genome (assembly: NCBI37/mm9). Word sizes 
of 6, 7, and 8 were analyzed with standard parameters. The maximal 
enrichment score (— logio p-value) observed for heptamers of each over- 
expressed miRNA was used for comparisons. For each dataset, the highest 
enrichment was observed in 3' UTR sequences within the threshold used to 
define "down-regulated" genes (—0.2 log 2 ratio), and the most represented 
words (7mers) exactly matched with the canonical seed (7mer-lA and 
7mer-m8) of the overexpressed miRNA (Fig. 6 B). Some of the miRNA- 
regulated gene lists (miR-22, miR-29, and let-Zb) also showed a weaker, 
but significant, enrichment of MREs in the CDSs (Fig. 6 C), which is consis- 
tent with previous evidence that a minor fraction of miRNA binding sites 
are located within the CDS (Lewis et al., 2005; Easow et al., 2007). None 
of the miRNAs showed enrichment in 5' UTR sequences. 

Ingenuity pathway functional analysis 

Significantly enriched functional classes (Figs. 1 C, 3 B, 7 B, and SI C) 
were identified through Ingenuity systems pathway analysis (Ingenuity Sys- 
tems), running a core analysis, and considering the enriched functions 
("Molecular and cellular functions" and "Physiological system development 
and function," P < 0.001) by using as input El A-regulated or miR-regulated 
genes. Each significant functional class corresponds to several signifi- 
cant (P < 0.001 ) gene ontology (GO) terms. Accordingly, for each class the 
minimum (MIN) and maximum (MAX) p-value of the corresponding GO 
terms was reported. The complete lists of all the significant functional 
classes for E 1 A-regulated or miR-regulated genes are reported in Table S2. 
To generate a network of cell cycle and DNA replication genes, the two 
most enriched functional categories of E 1 A Rb-dependent genes (cell cycle, 
206 genes, P = 2.60 x 10~ 57 ; and DNA replication, recombination, and 
repair, 2 1 4 genes, P = 3.65 x 1 0~ 33 ; Fig. 7 B) were selected. The com- 
bined (n = 270) list was used to build a network consisting of 1 38 proteins 
(nodes) and their interactions (edges) using the "connect" tool (with stan- 
dard parameters). The canonical pathways of "Gl/S checkpoint" and 
"Cell cycle control of chromosomal replication" were selected from the net- 
work and analyzed separately for mRNA-miRNA interactions. 

miRNA target enrichment analysis 

The mean number of either target genes or MREs (Fig. S5, A and B) for 
each pathway was determined as the ratio between the total number of ele- 
ments targeted (genes or MREs) within a pathway (N) and the number of 
miRNAs of each class (n), as follows: ALL, all the miRNA families included 
in the TargetScan 5.1 database (n = 373); PA, proliferation associated 
miRNA families (n = 13); DA, differentiation associated miRNA families 
(n = 26); and DA-7 miRNA families (n = 7). The following example illustrates 
the procedure: for all the 373 miRNA families listed in TargetScan, we 
found a total of 1,716 MREs (257 conserved and 1,459 nonconserved) 



in the "Gl/S checkpoint," and 780 MREs (42 conserved and 738 non- 
conserved) in the "Cell cycle control of chromosomal replication." Thus, 
on average, a single miRNA targets 4.6 Gl/S and 2.0 DNA replication 
MREs. Conversely, DA-7 miRNAs target a total of 57 Gl/S and 16 DNA 
replication MREs, with a mean of 8.14 and 2.3, respectively. Because all 
DA-7 miRNAs belong to the "broadly conserved" group of miRNAs, we 
repeated the analysis considering only broadly conserved miRNA families, 
as defined in TargetScan, (n = 86). We obtained almost identical results for 
all the analyses performed. To determine the significance of the difference 
observed, we generated 1 0,000 random lists (R-project software) of seven 
miRNAs (chosen from the 373 miRNA families listed in TargetScan), and 
we calculated how frequently a group of seven miRNAs randomly selected 
performed similarly (or better) than DA-7. 

miRNA-regulated genes overlap analysis 

Enrichment for Rb genes (Fig. 3 D) or differentiation genes (Fig. 4 B) among 
miRNA-regulated genes was performed according to Melton et al. (201 0). 
In brief, the enrichment was calculated by the ratio between the observed 
number of overlapping transcripts and the expected number of overlap- 
ping transcripts, as: [genes in overlap of miRNA-regulated group and Rb 
genes/all genes in miRNA-regulated group]/[all Rb genes/all genes used 
in analysis to generate miRNA-regulated groups]. As previously observed 
by Melton et al. (2010), regulated gene sets (either miRNA-regulated or 
target gene sets) are frequently enriched for highly expressed transcripts. 
Thus, to avoid a bias in the analysis, only highly expressed transcripts were 
used in the overlap (defined as those above the 50th percentile of average 
raw signal of all the Affymetrix expression profiles]. A similar pattern of 
results was obtained using all expressed transcripts, although all compari- 
sons appear more highly enriched. 

miRNA-mRNA interaction analysis on El A regulated genes 

We searched for DA-7 miRNA target genes (reported in Table S6) in the 
list of El A-regulated genes. To account for cotargeting of DA-7 miRNAs 
and to achieve enough statistical power, the analysis was performed using 
the union of all DA-7 miRNA targets. These targets were stratified accord- 
ing to the cumulative number of MREs for DA-7 miRNAs: 1 or more (3 1 3 
genes), 2 or more (173), 3 or more (101), and 4 or more (50). Then, we 
compared the regulation by El A of genes in the background list ("No 
target"] with that of genes targeted by DA-7 miRNAs (Fig. 9 B). We consid- 
ered the following gene classes: El A up-regulated genes, further stratified 
in Rb-dependent and Rb-independent genes. Statistical significance of the 
target distribution was calculated in JMP9 (SAS) with one-way analysis 
of variance (ANOVA; P < 0.001) and the Dunnett's t test (P < 0.05) as 
a post test, using the "No target" gene list as control group. Because the 
number of MREs usually correlates with the length of the 3' UTR region, 
the correlation of gene expression data with the number of MREs for DA-7 
was also adjusted for the length of 3' UTR regions by Multivariate Linear 
Model within JMP9. 

Clustering and statistical analyses 

Heat maps were generated by Java TreeView software for Mac OS X. 
Statistical analyses were performed within the statistical software JMP9. 
Significance of the differences was calculated using ANOVA (in the case 
of more than two groups] or Welch's t test. Categorical data were ana- 
lyzed with a contingency analysis within JMP (SAS]. P-values were calcu- 
lated with a Fisher's exact test. 

Online supplemental material 

Fig. SI and S2 show the El A-TD mRNA signature (GEX-MSC) and miRNA 
signature (MSC-miR) in MSCs, respectively. Fig. S3 shows that miRNA 
expression regulation during differentiation and upon El A expression 
occurs at a transcriptional level. Fig. S4 shows the dose response effects 
of individual miRNA expression on proliferation and the combinatorial 
effects of DA-7 miRNAs on protein levels of differentiation and prolif- 
eration markers. Fig. S5 shows the enrichment analysis of DA-7 miRNA 
targets within the Gl/S checkpoint pathway and the expression (as 
average copies per cell) of DA-7 miRNA families upon differentiation and 
E1A expression. Table SI contains the list of primers used in the lucifer- 
ase reporter assay. Table S2 contains details of the functional analysis 
performed on gene expression datasets. Table S3 contains details of the 
TRANSFAC analysis. Table S4 contains details of the El A-TD miRNA sig- 
nature. Table S5 contains details of miRNA- mRNA interactions within the 
Rb genes. Table S6 shows the list of miR targets. Table S7 contains details 
of the validation of miRNA-mRNA interactions for the Gl/S checkpoint 
pathway. Online supplemental material is available at http://www.jcb 
.org/cgi/content/full/jcb.201206033/DCl . 
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